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Abstract 

It is shown using experimental and numerical data that within the 
traditional inertial subrange defined by where the third order structure 
function is linear that the higher order structure function scaling expo- 
nents for longitudinal and transverse structure functions converge only 
over larger scales, r > rs, where rs has scaling intermediate between r\ 
and A as a function of R\. Below these scales, scaling exponents can- 
not be determined for any of the structure functions without resorting 
to procedures such as extended self-similarity (ESS). With ESS, differ- 
ent longitudinal and transverse higher order exponents are obtained that 
are consistent with earlier results. The relationship of these statistics to 
derivative and pressure statistics, to turbulent structures and to length 
scales is discussed. 

1 Introduction 

An important tool in understanding intermittency in turbulence has been the 
exponents £ p of power laws for the velocity structure functions. The longitudinal 
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structure functions are 

S%(r) =< (u(x + r) -u(x)) p >~ (1) 

where it and are in the same direction. Their measurement requires only a 
single hot wire probe that can, through the Taylor frozen turbulence assumption, 
determine one velocity component u as a function of the parallel spatial direction 
r and thus find the Sp at high Reynolds numbers. Using crossed-wire probes, 
one can also obtain reliable, high Reynolds number measurements of transverse 
structure functions (r), where it and ~f* are orthogonal, and their exponents 
£j '. There are also mixed structure functions containing both longitudinal and 
transverse components. 

The relationship between the £ p and intermittency jl) is in deviations of the 
the exponents £ p from their Gaussian or classical values of £ p = p/3, where 
£2 ~ 2/3 and C3 — 1 are expected for an energy cascade. In the presence of 
intermittency, £ p < p/3 for p > 3 is expected. Furthermore, it has generally 
been believed that all of the C^' T of a given order p should be the same in the 
infinite Reynolds number limit. This is closely related to the refined similarity 
hypothesis (RSH) which assumes that the only information that can affect the 
statistics on a given scale r is the fluctuations in the energy cascade e r through 
that scale. Details about either the large scale forcing or the dissipative struc- 
tures should be irrelevant in this picture. Measurements [g| have confirmed RSH 
as it relates dissipation to longitudinal statistics. 

More recently, moderate Reynolds numbers simulations and experiments 
[|| ^] have found, with the help of the extended self-similarity hypothesis 
that (J < f° r P > 3. This has been now been confirmed over all Reynolds 
numbers simulated || or observed to date. If longitudinal and transverse 
statistics are different, then the statistics of strain, that is dissipation e r , and the 
statistics of of vorticity, call it f2 r for its relationship to the enstrophy, should 
be different. Since £j < £^ , it implies that vorticity is more intermittent than 
strain. This possibility has been suggested by recent numerical results || that 
show that dissipation and enstrophy statistics scale as RSH predicts, but with 
separate distributions. This result implies that when the nonlinear coupling is 
written as a convolution in Fourier space that not only the velocity magnitude, 
but also the velocity phase, what makes the transfer of £l r different than e r , is 
important. Related to this, if the dissipation-dissipation correlation and its Q r 
counterpart go as 

< e x e x+r > - r" and < il x n x+r > - r^' (2) 

then different f2 r statistics implies that /z ^ [i. Experimentally, fj, « 0.25 
has been confirmed by the latest measurements and simulations |2^, |8| , 
although lower Reynolds number measurements tend to give /1 = 0.5. /j, ^ /x 
is found in one numerical result Ipfl . 
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It will be shown here that different longitudinal and transverse statistics 
suggest the existence of a crossover length scale within the inertial subrange. 
That is a length scale distinct from the large length scale L or a multiple of 
the dissipation or Kolmogorov scale 77, and maybe the order of the intermediate 
Taylor microscale A, where the Taylor microscale Reynolds number and A are 
defined as 

ux u 

R\= where A = - , 3 

v < (du/dxj 2 . > L ' Z 

R\ is the definition of the Reynolds number that appears to give uniform scaling 
in a variety of different flows and is related to the large scale Reynolds number 
by R\ ~ (R = UL/v) 1 / 2 . The Kolmogorov scale is related to A by 

n = (v/e) 1 /* ~ Ai?- 1/2 (4) 

where v is viscosity and e is the dissipation rate. 

In order to discern if there is some crossover length scale for some struc- 
ture function within the inertial subrange and determine if there are indeed 
different longitudinal and transverse exponents, a long inertial subrange and 
clean data are needed. Numerical simulations do not have this range, but if run 
sufficiently long give clean data and more flexibility. Observations can provide 
more dynamic range, but with more limited types of data. This paper will use 
both experimental and numerical data to try to present a more complete pic- 
ture than either measurements or simulations alone of what evidence there is 
for different longitudinal and transverse statistics and its implications. Exper- 
imental data will be used to indicate the high Reynolds number trends, then 
analysis from forced numerical turbulence in a box at R\ = 262 on a 512 3 mesh 
and R\ = 390 on a 1024 3 mesh will be presented, with most of the discussion 
related to the R\ = 262 calculation. It will be shown that if the same inertial 
range scaling analysis defined by the experiments is applied to the simulations, 
then the trends in the lower Reynolds number simulations are consistent with 
the experiments and clearly demonstrate a trend where simple scaling breaks 
down at many multiples of the Kolmogorov scale within the traditional inertial 
subrange, and that this scale appears to increase with Reynolds number. 

This paper will be organized as follows. First, constraints on and relation- 
ships between structure functions will be discussed. Recent results on structure 
function and pressure structure function scaling will be discussed in the context 
that different scalings for the same order imply the existence of a statistically 
signficant length scale within the inertial subrange where scaling properties of 
at least some structure functions could change. Next, there will be a discussion 
of existing experimental and observational results before showing new analysis 
of experimental data up to R\ = 3200. The analysis will use 5*^' T and S3 to 
define the minimum limits over which inertial subrange scaling analysis can be 
applied before the results for the ^ ' T for p = 4 and 6 are presented. Then the 
numerical analysis will be presented in a similar manner. 
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2 Constraints 



This section discusses kinematical constraints upon structure functions, related 
quantities, and what can be measured. There are only a few hard constraints 
for the scaling of small p structure functions. For p — 2, assuming isotropy and 
homogeniety and for a long enough inertial subrange, — Cj 1S expected. The 
classical value for £2 is 2/3 (||), which is closely related to assuming that the 
energy spectrum goes as 

E{k) = c k e 2/3 k- 5 ' 3 (5) 

where ck is the experimentally determined Kolmogorov constant and e is the 
equivalently either the rate of energy dissipation or the rate of energy transfer to 
small scales. With crossed wires, only one-dimensional longitudinal or transverse 
spectra can be found, but with a numerical simulation the full energy spectrum 
can be obtained in three-dimensions. An example from the 512 3 forced data set 
to be discussed is shown in Figure [l| This spectrum is found between a lower 
wavenumber, large-scale, forcing or integral scale L and a high wavenumber, 
small-scale, dissipation scale r\ = (v/e) 1 / 4 , with k v — 1/?]. 

The theoretical basis for predicting a -5/3 spectrum is the presumption of 
a local, uniform energy cascade. There is no constraint requiring -5/3 or the 
corresponding (2 = 2/3 for structure functions and corrections to the £2 = 
2/3 are claimed ||. However, no corrections to a -5/3 energy spectrum have 
been seen in numerical spectrum such as Figure [l] and for the largest range of 
scales observationally pj. How could two supposedly equivalent measures of 
turbulent scaling, the energy spectrum and the 2nd order structure functions, 
yield two inconsistent results concerning intcrmittcncy? A suggestion below is 
that dissipation range effects on S% are more persistent than thought. 

For p = 3, by balancing transfer and dissipation terms in the Karman- 
Howarth equation jy], |l2| , £3 = 1 is expected. This is the only fixed constraint 
on the £ p for isotropic, homogeneous turbulence regardless of intermittency and 
yields Kolmogorov's 4/5 law for 63 (r) 

S 3 (r) = ~er (6) 

However, linear in r behavior is never observed exactly (see discussion with 
figure ||). This effect has been quantified |l3| and it was shown that the peak of 
—Ss/r, after being compensated for the effects of forcing and dissipation, has 
a Reynolds number dependence like l s /L ~ i? A 3 ^°, whereas X/L ~ R^ 1 . This 
could suggest the existence of a dynamically significant position in the inertial 
subrange that does not scale with either the small Kolmogorov scale or the large 
scale and has been interpreted in terms of an enstrophy production argument 

While it is difficult to obtain a clean power law for p = 3, it is even harder 
to get convincing power laws when plotting S p versus r for p > 3. A device 
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Figure 1: Three-dimensional kinetic energy spectrum E(\k\) for the 512 3 forced 
calculation. The dashed line has an exponent of -5/3 and k v = fc max = 256. 
The dissipation spectrum e(fc) = vk 2 E{k) with a peak at k = 34 is shown. 
Resolution is determined by the ratio of the wavenumber of the peak of the 
dissipation spectrum (here k w k n /30 = 34) to fc max . Experience has shown 
that a ratio of at least 5 provides adequate resolution. 
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that has been used to determine higher order £ p is to make the assumption that 
even if — S3 ~ r is not exact, if it is assumed to be exact then much stronger 
power laws for the S p , p > 4, can be obtained by plotting S p versus —S3 rather 
than versus r. This is known as extended self-similarity or ESS ||. Another 
way of looking at ESS is that the S p are not true power laws all the way to the 
dissipation scale 77, but only for r > Aq, where A is large, and ESS allows one 
to extend scaling much closer to 77. However, there is a question of how much 
of this is a mirage because the ESS scaling exponents as r — > are trivially the 
classical values of p/3. The additional length scale Ar\ probably should have no 
dynamical significance since it is a multiple of 77, as will be discussed in section 

i 

For p > 4, there are no rigorous constraints that the longitudinal and trans- 
verse structure functions should have the same £ p . However, it has been shown 
[ p"5| that when all components of the velocity field it and all directions of d x 
are considered that there are four rotationally invariant combinations of fourth- 
order derivative correlations that can be written entirely in terms of the strain 
Bij and the vorticity and have been discussed in detail Q . 

F. = (15/7)<^> F eu = 3 <f e2 > (7) 
< e 2 > 2 < uj 2 X e 2 > 

< ujje i:i ej k LJ k > _ < a; 4 > 

f^eeu ~ 6 <L0 2 ><e 2 > *» ~ < ^2 >2 

Under the assumptions of homogeneity and isotropy F e is equivalent to the 
longitudinal fourth-order derivative flatness that can be measured with a 
single hot-wire probe and is related to the r — > limit of S%{r). To determine 
the other irrotational flatnesses as functions of R\ can only be done at low 
Reynolds numbers with complicated probes |fi"7f or simulations 16, [l8|], which 
suggest that each has its own scaling with Reynolds number. 

With a crossed-wire probe, higher Reynolds number observations of the scal- 
ing of flatnesses related to the r — > limit of the transverse and mixed structure 
functions could be determined. In addition to 



these are 



and 



F£ = Fuu =< uf A > I < u( A > 2 = F e (8) 



rpM _ p _ < > , Q x 

^4 — *utt — — t t (y) 

< u{ j >< U$ x > 



Fj = F tttt = 10 

< ul , > 2 



*2,1 

which can be related to combinations of the irrotational flatnesses (]?]). The scal- 
ing of F$f and Ff should be dominated by their most intermittent components, 
which the analysis of the irrotational components indicates is F^ [fL6| . 
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Generalized structure functions should involve not only more than one ve- 
locity direction, but also more than one spatial position and direction. All of 
the structure functions discussed so far involve different velocity components, 
but only one spatial position and direction. In this class would be included 
structure functions where the angle between 1? and d x is not or 90° fl9fl . 
The different rotationally invariant flatnesses (0) involve different velocity and 
spatial directions, but only one position, r = 0. It should be noted that a new 
type of structure function that uses two spatial positions in the same direction, 
has recently come into use in conjunction with fusion models [po[ . These are 
similar to the dissipation-dissipation correlation functions < e x e x+r > where 
the two distances would be (the derivative for dissipation) and the separation 
r between the two locations of dissipation e. 

While generalized structure functions could only be completely determined 
by numerical simulations or complicated probes at low Reynolds numbers, some 
insight might be gained by considering whether their number could be reduced. 
First, a way of systematically writing these down is needed p3|. However, 
unlike the way all fourth order derivative correlations can be expressed in terms 
of just four irrotational correlations, no simple reduction to a small number of 
fourth order structure functions has been found. What can be said is that if this 
number could be reduced, it would have to satisfy two rotational groups, one for 
position and the other for the velocity components [ p5[ . The full group would 
correspond to the spin plus angular momentum group in quantum mechanics. 
Even if such a reduction does not exist, one would expect some general properties 
to hold among all fourth-order structure functions. For example, if there is 
a subrange where the longitudinal and transverse structure functions, S^(r) 
and Sj(r), do have the same scaling, then maybe all fourth-order structure 
functions, including those related to fusion rules, are related in a similar manner. 
Or if there is a subrange where S^(r) and Sj(r), do not have the same scaling, 
then relationships such as the fusion rules should not apply. 

This leads us to three fourth-order structure functions measurable with 
crossed wire probes, each corresponding to one of the three derivative flatnesses 
(BHHJ). The longitudinal, mixed and transverse fourth-order structure functions 
are 

= Sou =< (u(x + r) - u(x)) 4 > 

Sf = S mt =< (u(x + r) - u{x)f{v{x + r) - v(x)) 2 > (11) 

Sj = S tttt =< (v{x + r) - v(x)f > 

whose general moment form is given by eqns. 13.83-13.84 of ^6|. However, the 
general moment form in no more fundamental than ( |Tl| ) in terms of the full 
rotational group. One would expect that if there are different scalings for the 
derivative flatnesses (0) under R\ , there should also be correspondingly different 
Cp' M ' T , with (p 1 and Cl more strongly dominated by the vorticity statistics. 
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3 Why two length scales? 



Another way of looking at the fourth-order velocity structure functions is to 
consider the second order pressure structure function f^M- This is related to 
a combination of fourth-order velocity structure functions ([ll]) by [^7) 

1 4 r°° 

P 2 (r) =< (p(x+r)-p(x)f >= --S^r) + -r 2 j y- 3 [^(y)+^(y)-65f (y)]dy 

( 12 ) 

+\f y- 1 [Sl(y)-^(y)]dy 

Assuming Sf (r) - Cf ,£f(r) - Co and 5 4 M (r) ~ (f then 

For the purposes here, the particular numerical prefactors are not important. 
The point to be made is that if the fourth-order longitudinal and transverse 
structure functions have different scalings, then P2 should have different scaling 
at the two ends of the inertial subrange. Assume that S% ~ r^ 4 and Sj ~ r^ 4 
with £4 > £j and the scaling for S^' 1 between these. Then for r small, P2 would 
be dominated by the scaling of 5j and for large r by S%, with some crossover 
length scale in the middle of the inertial subrange marking the separation be- 
tween two regimes for the scaling of P2. 

Related to P^if) is the pressure spectrum S p (k). It has been suggested Q 
that the dependence of Piik) on the spectral equivalents of the dissipation- 
dissipation correlation function < e x e x+r > and also < VL x VL x+r > and a cross 
correlation < £l x e x +r > places a kinematical constraint upon exponents for these 
correlations, and therefore the related fourth-order velocity structure functions. 
The relation used to show this is 

S P (k) = -^{E e (k) + E ( \k) E M (k)} (13) 

where E t (k), E il {k) 1 and E M (k) are the spectral versions of < e x e x + r >, v 2 < 
Q, x Cl x + r >, and v < Q, x e x + r > respectively. If each goes as 

E t (k) = C t e 2 k- 1 (kL)-^ (14) 

with corresponding f/ 1 and fi M for E n (k) and E M (k), it is then argued that 
unless fi e = /jP = fi M for all k there will be divergences in S p (k) as v — > 0. 
This is a reasonable conclusion for k — > 0. However, for k > k\ ~ 1/A where 
k*i/ 2 ~ 0(1), there is no need for this requirement. Therefore, there could be 
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a kincmatical constraint for < k < k\ requiring that fi e — fi = fj, x . Relating 
this to physical space, there would be a constraint that ( L — C, T — ( M as 
r — ► and R\ — > oo. This will be shown to be consistent with the analysis in 
sections (f|,[3|). For k\ < k < k v , there would be no constraint, so that it would 
be possible that /i c ^ ff 1 ^ /i M and C L / C T 7^ C M as r — * which is also 
indicated by the analysis here. What would be most satisfying would be if the 
pressure spectrum itself showed a clear break near k = 1/A. This has now been 
found in the pressure analysis of the 1024 3 data to be discussed here Q . 



4 New experimental evidence 

There have been numerous experimental and observation studies of turbulence 
designed to determine the scaling of structure functions and related measures 
of intcrmittency such as the dissipation-dissipation correlation function (|J). 
However, careful examination raises several questions. First, while the high- 
est Reynolds number observations do find that C T (p) 7^ C L (p); when ESS is 
used, compared with lower Reynolds number results, the difference is noticeably 
less. This suggests that in the very high Reynolds number limit the difference 
could go to zero. On the other hand, ESS might be giving a false impression of 
good scaling since the ratios Sp(r)/Sp(r) as functions of r do not show good 
scaling behavior to as small a scale as plotting with ESS does. A recurrent 
limitation is that scales below A in experimental and observational analysis are 
usually not shown, perhaps because it is felt that the small scales are less reli- 
able. For example, the dissipation-dissipation correlation function exponent is 
expected to be /i « 0.25 from the observed value of ^ pa 1.78. This is con- 
firmed by several high Reynolds number measurements (p0[ p2[ and references 
therein), but only for r > A. It is noted |22j that lower R\ experiments tend 
to give (i Pa 0.5. One moderate Reynolds number (Rx = 400) experiment for a 



circular jet 23 that shows \i pa 0.2 for r > A, also shows a breakdown in scaling 
of < e x e x+r > for r < X that could be with the origin of /x s=s 0.5 at low R\. 

The new experimental analysis comes from data from two experiments at 
the Moscow wind tunnel that has been used to investigate a number of funda- 
mental issues involving turbulent spectra and structure functions Jl(], . The 
two experiments are for a mixing layer ML at i?A=2100 and a return chan- 
nel RC at i?A=3200, which are the highest Reynolds number laboratory data 
sets available with both longitudinal and transverse velocities. The advantage 
of an experiment over atmospheric observations is that an experiment offers 
more controlled conditions, which could be especially important for determin- 
ing transverse structure functions since any inherent anisotropy in the flow could 
affect their values. With respect to this, an important point to remember in the 
following discussion is that the ML data set is very anisotropic and RC data set 
is nearly completely isotropic. Therefore, only the RC, R\ = 3200 experimental 
data can be directly compared to the isotropic numerical data to be presented. 
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Figure 2: The normalized third-order velocity structure functions. 4> u $ = — < 
Auf > /(< e > r). x , ML; o, RC. The vertical down and up arrows indicate the 
classical inertial-range bounds (r v and tl) for ML and RC, respectively. The 
isotropic value of 4> u ,3{r) — 0.8 is shown by the solid line. 
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The ML data set is included to show that most of the relevant properites also 
appear when there is anisotropy. Extended self-similarity will not be used in 
order to emphasize the regimes with simple scaling and where this breaks down 
within the inertial subrange. It will be shown that the scalings of all of the 
measured higher-order structure functions change within the inertial subrange. 
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Figure 3: The normalized second-order velocity structure functions. 4> s 2 =< 
As 2 . > /(< e > r) 2 / 3 . o, s = u, <p Ut2 = 4>2 I x , s = w, <t> w $ = 4% ■ V = 0.00018 
and 0.00014 for ML and RC respectively. Arrows are r v and r^. The vertical 
arrows indicate the classical inertial-range bounds. 

Figure ^| shows the third-order longitudinal structure functions — Ss(r) di- 
vided by r for ML and RC. Excluding using the absolute value, all odd order 
transverse-only structure functions are zero. Figure ^ uses arrows to indicate 
the inertial subrange based upon the regime over which — S3 (r)/r is constant for 
ML and RC. For ML and RC respectively, at the large scale arrows are lengths 
tl = L/5 and L/6 and at the small scales lengths r v — 307/ and 25r/. Fig- 
ure H shows the second-order longitudinal S^ir) and transverse 5*J(r) structure 
functions for ML and RC. For S^r), at large r there is scaling out to an Tl con- 
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sistent with tl from S3 (r) in Figure g. In addition, there is a new large length 
that we will call r 4 / 3 < tl over which the isotropic relation Sjir) rj (4/3)5f(r) 
holds. r 4 / 3 is smaller than the range over which the longitudinal second order 
structure function S% has good scaling behavior, which is roughly out to Tl, 
but is consistent with the maximum r for which the transverse second order 
structure function S% has good scaling. For ML and RC, r 4 / 3 = L/30 and L/25 
respectively. The small length scale r v is roughly where the one must begin 
to apply extended self-similarity, ESS, (section ||) if one is to get good scaling 
relationships for the higher order longitudinal structure functions down to the 
Kolmogorov scale r\. The reason r n is being introduced is to clearly indicate 
that any new length scale between r v and r 4 / 3 or tl is fully within the incrtial 
subrange. While r 4 / 3 < tl at large scales, the range over which the 4/3 rule 
for S 2 ' (r) seems to fit at small scales extends to r < r v . This would be con- 
sistent with how extended self-similarity extends scaling regimes more into the 
dissipation regime. 

Figures |](a-d) plot the logarithmic derivatives of the fourth and sixth-order 

longitudinal and transverse structure functions. For fourth and higher order 

structure functions there are no isotropy relationships that would require that 

Sp have the same scaling as over the entire inertial subrange, only the 

pressure spectum argument (section ||) that Cj should equal £4 for r > A. 

Consistent with the pressure spectrum argument, there is a regime in Figures 

[|(a-d) where the higher order structure functions do have the same scaling, 

between two lengths, r 4 / 3 and a new small length scale that we will call r$ > r n . 

Referring to Table 1 for ML (R x « 2100) and RC (R x « 3200), r s = 60f? and 

1/2 

50?7 respectively while A = R x rj rj 46/y and 56rj respectively. So rs ~ A, but 
it is also roughly twice the value of r, r For r < rs, both £j and increase 
rapidly from their r > rs constant values, with ^ increasing the fastest. 

This analysis shows a regime within even the rather strict definition of a 
measurable inertial subrange defined by r v < r < r 4 / 3 where universal scaling 
of longitudinal and transverse structure functions is found and therefore as- 
sumptions of refined self-similarity might apply. In addition, there is a regime 
for r < rs where both the longitudinal and transverse scaling exponents diverge 
from constant behavior. Praskovsky's analysis shows that if scaling functions 
were fitted over the entire inertial subrange, that different exponents for the 
longitudinal and transverse structure functions would be found that would be 
consistent with recent experiments and simulations [|[ m. 

5 New numerical evidence 

To determine more clearly how significant the differences between the longitu- 
dinal and transverse structure functions are, it is necessary to see some trends. 
For this purpose we now move to analysis of a 512 3 forced calculation of isotropic 
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Figure 4: The local slopes Cp (r) = ^log^' 2 (r)]/dlog(r), p = 4, 6. s 



for Cp and s = w for Cp • °, C«,p0) = Cp ( r ); x > C«-,pM = Cp M; *, R P (r) 
?p L (r)/Sj 



Sp(r)/Sp(r). The outer small vertical arrows indicate the classical inertial 



range bounds, r n and r^, while the inner large vertical arrows indicated the 
new bounds, rs and r 4 / 3 , where Cj ~ Cp ■ 



13 



turbulence in a periodic box with R\ = 262 and another 1024 3 calculation with 
R\ = 390. The cleaner numerical data also allows one to more directly compare 
S% and S2 and to apply extended self-similarity. 

The simulations are classic turbulence simulations in a periodic box with the 
lowest band of wavenumber modes forced to have constant energy by Gaussian 
white noise. Statistics were taken several eddy turnover times after any large 
excursions in the dissipation due to the initial conditions have dissappeared and 
once the dissipation rate had settled to the point where it was varying by less 
than 5% on the timescale of several eddy turnover times. Experience has shown 
that this type of Gaussian white noise forcing yields stable statistics over a single 
eddy turnover time. The 512 3 simulations were not dealiased and the statistics 
represent an average over 40 large scale eddy turnover times done on the Cray 
T3D of the IDRIS Institute in Orsay, France. The 1024 3 calculations were done 
on a Fujitsu V005999/56 at the Nagoya University Compution Center and 
represent an average over 1.1 large scale eddy turnover times. How well the 
1024 3 calculation satisfies isotropy relations for S% ' T (r) and Ss(r) is discussed 
elsewhere |5l| . Figure [j] is the three-dimensional kinetic energy spectrum for 
the 512 3 calculation, showing a clear -5/3 regime. The Kolmogorov scale 77 is 
half the mesh size r\ = Ax/2. This is considered adequate resolution because 
the peak of the dissipation spectrum, roughly where the -5/3 regime rolls over 
into the dissipation regime, is at k w k n /30 where k n = 2it /r). 

In the experiments above, the original definition of an inertial subrange 
placed the smallest scale useful for analysis at about r v — 25r/, roughly the 
scale associated with the peak of the dissipation spectrum. This definition of 
r n would make it greater than A for both simulations, which for R\ = 262 is 
A w I677 and for R\ = 390 is A w 20ry. Therefore, a quantitative means of 
choosing r v must be determined that can be applied equally to simulations at 
lower Reynolds number and the R\ = 3200 experiments. A difficulty in choosing 
a method is that the simulations are not at high enough a Reynolds number to 
exhibit either a clearly defined linear regime in S3 or a long regime where the 
isotropy relationship between 5^ and 5j is obeyed. 

The problem is illustrated by Figure ||, which plots — S3(r)/(.8er) and the 
normalized form with the viscous correction 

, N -S 3 (r) + Qv^l , x 

D LLL (r)= 3W ger 9r ■ (15) 

While the span over which Dlll{t) is 1 indicates that there is an energy cascade, 
an inertial subrange over which — S^r)/ .8r is flat as in the experiments in Figure 
H is not seen. Therefore, arbrtrarily a inertial subrange will be chosen as those 
r where — Ss(r) / (.8er) > 0.8. This is the consistent with how r v was chosen 
for the RC experiment and gives consistent values of r v for the 512 3 and 1024 3 
simulations and the RC experiment. Based upon Figure ||, an inertial subrange 
can be defined between r/r) =19 and 193. 
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Figure 5: -S 3 (r)/(.8er). Also shown is {-S 3 (r) + 6vdStf(r)/dr)/(.8er). An 
arbitrary constant line is drawn to indicate where an inertial subrange is defined, 
r is with respect to the size of the box (2tt) 2 
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Let us now use S%{r) and S'J(r) to consider this definition of an inertial 
subrange using the isotropic relationship in 3D between the transverse and lon- 
gitudinal structure functions. 

Sl(r)^S^(r) + r -dS^r)/dr (16) 

In the inertial subrange where 62 (r) ~ r 2 / 3 , from jl^ ) one gets 5j(r) ~ 
(4/3)S%(r). As r -> 0, S 2 (r) ~ r 2 , so one gets Sf (r) - 2Sf(r). Figure | 
shows the second-order longitudinal and transverse structure functions divided 
by r 2 / 3 . S-f (r)/r 2 / 3 is plotted twice, first divided by 4/3 to demonstrate how 
well the inertial range relation is obeyed and then divided by 2 to demonstrate 
the approach to the dissipation range. 

The 4/3 rule in Figure || is only approached at the largest scales, for r 
greater than where — Ss(r) / (.8er) is greatest in Figure^, which could be called 
''max- For r < '"max m Figure H the gradient of S-f (r) is slightly steeper than 
the r 2 / 3 prediction, which would be consistent with measurements of a small 
correction over many years that is usually interpreted as due to intermittency. 
However, 5*J(r) does not show this correction, and since it is S^ir) which forms 
the major portion of the energy (4/5ths), it should not be surprising that the 
energy spectrum (Figure |l|) is very close to -5/3. If the slope in S% (r)/r 2 / 3 is a 
dissipation range effect, then one could compensate for this in analysis of higher 
order structure functions by using a new ESS variable based upon 8% (r) 

R$ = (S^(r)f/ 2 /e, (17) 

rather than — 5a(r). This is done below. 

For direct comparison with the experiments, Figures 0a and (?]b show the 
logarithmic derivative of the 4th and 6th order structure functions for R\ = 262 
against r and Figure || shows S^' M ' T (r) for R\ — 390. As in Figure Q there 
is a span we will define as tl > r > r$ where the slopes of the longitudinal 
and transverse structure functions are constant, but in this case they are not 
identical. In the experiments, r$ was chosen to be the first r (from below) 
where Cj( r ) Cp( r )- Since there is no regime where Cj( r ) — (pi r ) m the 
simulations, the choice of r$ is more subjective for the simulations. We have 
chosen r$ to be where a line of constant (4 = 1.28 at large r would meet a line 
with a logarithmic dependence through <^J at small r. The results are given in 
Table 1 and are similar to the experiments in that r$ is greater than r^, the 
lower limit that was defined for inertial range behavior in S3. 

What is particularly similar to the experiments in Figure [?] is that for r < r$, 
all (p ,T increase and this increase is greater for than for £j. There is also a 
clear trend where the difference between (J and (p defined as 

5 p = min{C p L - Cj} (18) 
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Figure 6: S% (r)/r 2 / 3 (line) and 5'J(r)/r 2 / 3 plotted twice. Once divided by 
4/3 to show the inertial subrange (dash) and once divided by 2 to show the 
dissipation range (dot-dash). 
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is decreasing as R\ increases, pointing to the experiments where Cj ~ f° r 
r > r s . 

Applying ESS with R2 brings the slopes back down for r < rs, long regimes 
for scaling S^(r 3 ) and Sg (r 3 ) appear, and the differences between £j and £p 
appear more clearly. For the R\ = 262 simulation, Figure ^| shows different 
C,^ and £p that are consistent with earlier work in the same Reynolds number 
regime jij, |(|. However, the long regime of nearly constant should now 
be considered an artifact of using ESS and the need to apply ESS seems to be 
intimately tied to the new length r$, where rs is not simply a constant times 
eta. 
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Figure 9: and £j for R\ = 262. Values are taken at r/77 = 188. SL is the 
She-Leveque formula (32), which is an excellent fit to the observations, p/3 is 
the classical prediction. 



6 Discussion 

In Kerr ]lq| , in addition to velocity statistics, equivalent passive scalar and 
mixed velocity-scalar statistics were calculated. It was noted that the scaling of 
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the derivative flatness for a passive scalar 9 

F 99 a = F vg = (9/5) < (Vfl) 4 > /(< (V0) 2 >) 2 (19) 

was similar to the vorticity flatness F u , that is z\ |33) and z% were both much 
larger than the exponent z\ for — F e . There was also a strong anti- 
correlation between the scalar derivative and vorticity, that is 

Fve^ =< ($0 ■ rtf > I < (V0) 2 >< oj 2 > < 1 (20) 

Taking this analogy between the statistics of the scalar gradient and vorticity 
a step further, if the scaling of Fve could be used as a tool for determining the 
scaling of F u , then high Reynolds experiments for temperature statistics Q 
could have been implying greater scaling exponents in the transverse derivative 
correlations and greater deviations from classical structure function exponents 
long before the new work with crossed-wire probes. However, without more 
theoretical understanding and corroborating evidence from velocity structure 
functions, these hints were not studied further. 

The theoretical, experimental and numerical discussion here replaces that 
speculation with moderate to high Reynolds numbers experimental and numer- 
ical data about the behavior of the transverse structure functions. The impor- 
tant points are that for p > 3 that £j w C£ for r > rs as R\ — ► oo, where 
rs is further into the inertial subrange than expected, and that none of the £ p 's 
show simple scaling behavior for r < rs- Table 1 shows $4, <5 6 , r„, rs and A 
for the numerical cases of R\ = 262 and R\ = 390 and the experimental cases 
ML (R x = 2100) and RC (R x = 3200). Between ML and RC, mixing layer 
and return channel, the trend for r^ and rs is opposite that between the much 
lower Reynolds number simulations and the experiments. This is probably just 
a reflection of the influence of anisotropy for ML. The overall trend supports 
the existence of a length scale rs that is much larger than 77, with the scaling 
of rs as a function of R x intermediate between that for rj and that for A. 

The theoretical discussion showed that this second small length scale would 
not be expected if the energy cascade and the refined similarity hypothesis 
were controlled only by the statistics of the dissipation e. For there to be a 
dynamically significant length scale within the inertial subrange, there must be 
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something in addition to e controlling the cascade, the fundamental dissipation 
mechanism must involve two length scales, or both. This paper addresses only 
one kind of higher order statistic, the structure functions. Clearly a thorough 
analysis of all higher order statistics, including pressure and the dissipation- 
dissipation correlation function, needs to be done on available measurements 
and simulations from the point of view of determining whether some length scale 
of the order of the Taylor microscale has a role for them also. This has been 
done for pressure spectra for the 1024 3 data set used here |29| and the results are 
consistent with the existence of such a length scale separating spectral regimes 
of -7/3 and -5/3. 

Either numerically or observationally, consistent conditions over a wide range 
of Reynolds numbers are necessary if any conclusions are to be drawn. This is 
difficult to obtain with atmospheric measurements. As an example, for one 
series of atmospheric measurements over a wide range of Reynolds number, 
< e x e x -i- r > has been determined over the entire inertial subrange JlO| . For 
larger r, there is convergence to /i w 0.2, but there is an enormous scatter 
between different measurements at smaller r so that any R\ dependence in the 
equivalent of rs for < e x e x+r > would be difficult to determine. Therefore, well- 
controlled high Reynolds number experiments would be very useful. This could 
also provide a motivation for doing a careful 2048 3 forced simulation, which is 
now feasible. 

If confirmed, what could be the dynamical significance of this crossover 
length scale beyond just being an average between the integral scale and the 
Kolmogorov scale? To date, no importance has been attached to the Taylor 
microscale A beyond that of an average. When the first visualizations of vor- 
tex filaments were done |L6]], one way of characterizing them was they had a 
width the order of the Kolmogorov microscale and a length the order of the 
Taylor microscale. However, in a 64 3 or 128 3 DNS, it would be impossible to 
determine whether the length was A, a fraction of the size of the box, or just a 
multiple of r\. That is, the the radius of curvature might just be a multiple of r\. 
The highest resolution visualizations of isotropic, homogeneous turbulence [p5| 
would support this scenario. That is, vortex filaments are observed, but they 
do not snake through the entire domain and instead have a length that appears 
to be a only a multiple of r\. However, it can be argued that these are hy- 
perviscous calculations that are predisposed to shortening the vortex filaments. 
Furthermore, statistical models based strictly upon vortex filaments |32, Q do 
not seem capable of producing different longitudinal and transverse scaling in 
the high Reynolds number limit. 

Therefore, it seems that some other type of dynamical object besides sim- 
ple filaments would be needed if a theoretical basis for a second dynamically 
significant length scale is to be given. The only dynamical structure that has 
been identified in cither turbulence or idealized calculations of Navier-Stokes 
and Euler that self-generates two small length scales is the structure found in 
the interaction of two anti-parallel vortex tubes as the peak vorticity appears 
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to be developing a singularity |57|. However, the spectrum of this structure is 
fc~ 3 , nowhere near /c~ 5 / 3 . So until high Reynolds number, very highly resolved 
calculations are done for Navier-Stoke vortex reconnection, any connection be- 
tween the properties of this structure and the scaling properties of turbulence 
is pure conjecture. 
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